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Introduction 


High  speed  digital  computers  have  had  a  significant  impact  on  the  rough  surface 
scattering  community  because  they  provide  a  new  and  important  problem  solving  tool. 
As  an  example,  computer  simulations  have  led  to  a  much  better  understanding  of  the 
root  causes  of  the  phenomenon  known  as  enhanced  backscattering  from  rough  surfaces 
[1].  However,  to  date,  it  has  only  been  possible  to  use  the  full  power  of  the  computer  on 
one  dimensionally  rough  surfaces  because  the  two  dimensional  rough  surface  scattering 
problem  is  still  too  large  to  be  so  solved.  There  is  work  being  done  with  techniques 
which  might  make  the  two  dimensional  problem  more  tractable  on  a  computer  [2]; 
however,  such  improvements  remain  to  be  thoroughly  demonstrated. 

This  state  of  affairs  leads  to  the  rather  obvious  conclusion  that  there  is  still  a 
need  for  further  analytical  work  to  reduce  the  numerical  difficulty  of  the  problem  to  a 
more  tractable  level.  The  intent  here  is  to  isolate,  separate  out,  and  perform  as  much  of 
the  solution  as  possible  using  analytical  techniques;  numerics  will  only  be  required  for 
those  parts  of  the  problem  where  analytical  approximations  are  not  possible  or 
jeopardize  the  accuracy  of  the  solution.  The  key  ingredient  in  such  a  solution  is  the 
development  of  an  approach  which  is  very  robust  in  its  capabilities.  For  example,  one 
would  not  want  to  start  out  with  a  boundary  perturbation  type  solution  in  which  both 
the  surface  height  and  the  surface  slopes  are  required  to  be  small  if  the  solution  were 
going  to  be  applied  to  surfaces  with  large  slopes.  Two  approaches  have  recently  been 
introduced  and  are  being  developed  which  appear  to  hold  the  promise  of  significantly 
improving  our  capabilities  with  regard  to  predicting  the  scattering  from  two 
dimensional  rough  sur  aces.  The  first  of  these  is  the  work  of  Rodriguez  and  Kim  using 
what  they  call  the  un  fied  perturbation  expansion  approach  [3];  this  technique  is  based 
on  a  small  momentum  transfer  approximation  developed  earlier  by  Rodriguez  [4].  The 
other  approach  is  based  on  the  method  of  smoothing  [5,  6,  7],  a  demonstration  of  its 
ability  to  deal  with  surfaces  having  arbitrary  surface  height  derivatives  for  small  heights 
[8],  and  a  possible  extension  to  surfaces  having  larger  heights  [9].  It  should  be  noted 
that  both  of  these  approaches  hold  the  promise  of  extending  our  prediction  capabilities 
well  beyond  those  of  existing  asymptotic  solutions.  The  remainder  of  this  report  will  be 
concerned  with  the  solutions  which  are  based  on  the  method  of  smoothing. 

Whenever  new  approximate  solutions  are  obtained  for  complex  problems  such  as 
scattering  by  randomly  rough  surfaces,  it  is  wise  to  compare  these  with  existing 
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approximations.  This  is  clone  for  two  reason.  First,  it  is  essential  to  compare  with 
known  benchmarks  in  order  to  check  the  accuracy  of  the  new  solution  and  to  see  where 
it  provides  improvement.  A  second  purpose,  which  is  frequently  overlooked,  is  to 
generate  new  ideas  as  to  how  the  new  approximation  may  be  improved  even  further.  It 
is  this  latter  reason  that  forms  the  basis  for  this  report.  The  basic  intent  is  to  compare 
what  has  been  called  the  normalized  first  order  smoothing  (NFOS)  approximation  with 
the  conventional  composite  surface  scattering  model  to  (1)  first  determine  what 
improvements  the  latter  suggests  for  the  former,  and  (2)  to  then  demonstrate  just  how 
important  these  improvements  are.  The  final  result  will  be  a  new  approximate 
scattering  result  which  essentially  extends  the  existing  composite  surface  model  to  a 
new  class  of  surfaces.  Unlike  the  conventional  composite  surface  model,  this  new  result 
will  also  contain  a  well  defined  methodology  for  extending  its  capabilities  in  the  sense  of 
making  it  apply  to  more  complicated  surfaces.  This  is  a  feature  that  the  conventional 
composite  surface  model  never  exhibited. 

Before  presenting  the  main  topic  of  this  report,  it  is  first  necessary  to  review  the 
first  order  smoothing  (FOS)  approximation  and  its  capabilities  and  limitations.  Of 
particular  note  is  its  ability  to  reduce  to  the  Rice  boundary  perturbation  approximation 
provided  the  surface  height  is  very  small  in  terms  of  the  electromagnetic  wavelength 
and  the  angle  of  incidence  is  not  too  close  to  grazing.  Of  equal  importance  is  to  show 
why  FOS  begins  to  breakdown  as  the  surface  height  approaches  or  exceeds  the 
electromagnetic  wavelength.  Identification  of  the  failure  of  FOS  as  the  height  becomes 
large  suggest  a  normalization  scheme  to  overcome  this  limitation  and  this  leads  to  the 
normalized  first  order  smoothing  (NFOS)  approximation.  By  analyzing  the  class  of 
surfaces  to  which  NFOS  should  apply  and  comparing  these  with  the  surfaces  to  which 
the  conventional  composite  surface  model  applies,  a  means  for  improving  the  NFOS 
approximation  even  further  is  developed.  It  is  this  improved  NFOS  model  which 
extends  the  basic  smoothing  technique  to  a  class  of  surface  not  heretofore  addressable. 
The  final  result  is  a  smoothing  based  scattering  approximation  that  is  applicable  to  the 
type  of  multiscale  roughness  surfaces  that  are  common  in  nature. 

Background 

The  specific  problem  to  be  considered  is  that  of  scattering  of  an  incident  plane 
electromagnetic  wave  by  an  arbitrarily  roughened  perfect  electric  conductor  occupying 
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the  half  space  defined  by  z  <  C(x,  y).  The  roughness  height  is  taken  to  be  zero  mean 
about  the  z=0  plane,  i.e  ,  (C)  =  0,  and  is  assumed  to  be  a  statistically  homogeneous 
process.  This  latter  constraint  implies  that  the  surface  statistics  are  independent  of  the 
point  (x,y)  where  they  are  being  measured.  The  region  above  the  roughened  conducting 
half  space,  e.g.,  z>C  (x,y),  is  taken  to  be  free  space.  The  incident  plane  wave  is 
traveling  in  the  direction  specified  by  the  unit  vector  k,-  and  ks  indicates  the  direction  of 
scattering,  see  Figure  1. 

The  current  induced  on  the  conducting  surface  Js(r)  obeys  the  magnetic  field 
integral  equation,  e.g., 

J.(f)  =  J,''(f)  +  2fi(f)«  Arx  Js(f0)dS„  (1) 

which,  after  multiplying  by 

\A  +  (  v  02  exp  (jk„  C), 


becomes 

J(0  =  J*(f )  +  Lr(r,  r0)  •  J(r0)  (2) 


where 


VC  = 


dC  ~  ,  (K  ~ 
d^x  +  dF 


+  z 


(3a) 


ksjr  —  k0  ks  •  z 


(3b) 


“(?)  =  (^5f-^y  +  z)/V/l  +  (V02  =  \/l  +  (VC)2  N(r)  (3c) 

g(|Arj)  =  exp(-jk0|Ar|)/47r|Ar|  (3d) 

Ar  =  (f-f0)/|r-fo|  (3e) 


and  kQ  =  2n/X0  and 


.J(f)  =  Js(f)  \/l  +  (  VC)2  exp(jkS2C) 


(4a) 
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Figure  1 .  Scattering  geometry 
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Jl(f)  =  2N(r)  x  H*(r)  exp(jkS7A()  (4b) 

oo  oo 

L  =  1 1  dx0dy0  =  J  |  dp0  (4c) 

-  oo  —  oo 

r  =  -2  exp  (jks*A<)  N(r)  x  (Ar x  (4d) 

AC  =  C  (p)-Upo)  (4e) 


In  the  above,  n(r )  is  the  unit  normal  to  the  surface  at  the  point  r  =  p  +  Cz  on  the 
surface,  d(/dx  and  d(/dy  are  the  x  and  y  surface  slopes,  g  is  the  free  space  Green’s 
function  evaluated  on  the  surface,  and  H*  is  the  incident  magnetic  field.  It  should  also 
be  noted  that  r(r,f0)  is  a  3x3  matrix  operator.  The  reason  behind  manipulating  (1) 
into  (2)  is  that  a  simple  Fourier  transform  of  the  currents  appearing  in  (2)  essentially 
yields  the  far  zone  fields  scattered  by  the  rough  surface.  In  fact,  previous  papers  dealing 
with  the  smoothing  method  [8,  9]  have  found  it  advantageous  to  deal  with  the  transform 
of  (2);  however,  in  this  paper  attention  will  be  directed  toward  the  current  equation. 
The  reasons  for  this  will  become  obvious  later. 

The  basic  idea  behind  smoothing  is  relatively  straightforward.  First,  the  current 
is  split  into  a  zero  mean  and  a  fluctuating  parts,  i.e. 

J  —  (J)  +  6  J  (5) 

Quite  apart  from  the  approximations  this  leads  to,  (5)  is  a  natural  decomposition  in 
that  (J)  gives  rise  to  the  specular  or  coherent  scattered  field  while  63  leads  to  the  diffuse 
or  incoherent  scattered  power.  Substituting  (5)  into  (2)  and  successively  taking  the 
average  and  fluctuating  part  of  the  resulting  equation  generates  two  coupled  integral 
equations  for  (J)  and  <5.1.  For  example,  the  equation  for  63  is 

63  =  [1  -  L(r)] _  1  (5r.(J)  +  £[L6r.6J]}  (6) 

In  (6),  6V  is  the  fluctuating  part  of  the  propagator  in  (2),  and  (1  -  L(r)]  “  1  is  an  operator 
m  coordinate  space  but  becomes  a  multiplicative  function  in  transform  space.  The 
specularity  of  the  average  scattered  field,  which  has  been  demonstrated  quite 
independent  of  the  smoothing  method  [10],  implies  the  following  simplification 
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<J)=j  exp 


(?) 


where  =  kj  -  kiz  £  and  j  is  a  vector  number  which  is  independent  of  any  of  the 
surface  coordinates.  Substituting  (7)  into  (6)  yields 

63  =  [1  -  L(r>]  "  1  { <5  J  *  +  Sf(r,  -  ki()  .f 

-f  6[L6T«  6J]>  (8) 

where  the  tilde  over  6T  denotes  the  Fourier  transform  from  p0  space  to  klt  space. 

For  most  surfaces  of  interest,  a  good  approximation  for  j  is  j  *  which  is  defined  by 

jl  =  2H0  (Nxh,-  exp[j(k„-klz)c])  (9) 

where  H0  is  the  amplitude  of  the  incident  magnetic  field  and  hi  is  its  vector  direction. 
For  Gaussian  surfaces,  the  height  and  slopes  at  a  common  point  on  the  surface  are 
uncorrelated  so  (9)  becomes 

j  ‘  =  2H0  (Nxh,-)  (exp[j(ksz  -  k, -,)<})  (10) 

The  average  of  the  exponential  height  term  is  recognized  to  be  the  characteristic 
function  of  the  surface  height.  Substituting  (10)  in  (8),  ignoring  the  operator  [1  -L  (r)] 
since  it  has  already  been  assumed  to  be  small  (in  replacing  j  by  j1),  replacing  5r  by  r, 
and  retaining  only  the  new  Born  term  in  the  resulting  integral  equation  yields  the  first 
order  smoothing  (FOS)  approximation; 

*J(1W'  +  2H0  (exp[j(kM-kizK])  r  .(Nxh1)  (11) 

It  can  be  shown  that  6j*  is  the  term  that  is  most  important  about  the  specular 
direction,  i.e. 

E'|,L-Tr7-rs  (12) 

!  kjt  i 
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while  the  second  term  is  the  major  contributor  far  away  from  the  specular  direction.  In 
addition,  it  can  be  shown  (although  not  easily)  that  (11)  leads  to  the  classic  Rice 
boundary  perturbation  result  provided  that  k0C<<l  and  the  incidence  angle  is  not  too 
close  to  glazing  incidence.  The  first  restriction  permits  ignoring  the  height  dependence 
in  6r  while  the  second  insures  that  the  terms  which  are  ignored  are  indeed  small  relative 
to  the  ones  retained  in  the  simplification  process.  It  is  interesting  to  point  out  that 
simplifying  (11)  to  the  Rice  limit  does  not  require  small  slopes:  this  would  tend  to 
indicate  that  the  Rice  limit  is  probably  more  robust  than  classical  derivations  would 
tend  to  indicate.  Certainly,  there  are  cases  in  the  literature  where  the  Rice  limit 
appears  to  “work”  when  the  classical  restrictions  would  suggest  otherwise  [11]. 

Interestingly  enough,  the  FOS  appr  ximation  in  (11)  clearly  show's  why  it  fails  as 
the  surface  height  becomes  large  relative  to  the  EM  wavelength.  More  specifically,  as 
the  surface  roughness  height  becomes  large  t  ‘  e  characteristic  function  for  the  height  in 
(11)  decreases;  for  a  Gaussian  surface,  the  decrease  is  most  rapid.  This  behavior  is 
clearly  incorrect  because  the  increase  in  height  can  be  brought  about  by  the  addition  of 
a  randomly  elevated  planar  component  to  the  surface  roughness  which  has  no  horizontal 
structure  and  yet  it  leads  to  an  attenuation  of  the  wide  angle  scattering  contribution  in 
(11).  The  difficulty  here,  as  identified  and  remedied  in  [9],  is  that  smoothing  is  only 
applicable  to  surfaces  having  a  relatively  small  height.  The  remedy  for  this  situation  is 
to  return  to  the  original  current  integral  equation,  (2),  and  normalize  this  equation  by 
the  height  dependent  phase  factor  appearing  in  the  term  J*  in  (2).  This  normalization 
leads  to 


J„  =  J  exp[-j(ksz  -  ktr)c] 

(13.i) 

■V  =  2NxH*  exp[jkj-,C] 

(13b) 

r„  =  r(f,  r0;  k«  -  k-J 

(13c) 

hi  =  -hi*  +  L  !\i  * 

(14) 
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Applying  smoothing  to  this  normalized  equation  yields  the  following  integral  equation 
for  the  fluctuating  normalized  current; 

*J«  =  [l  -  L<r„)] " 1  +  srn .  fn  +  6[L6rn .  *J„]}  (15) 

Assuming  that  jn*  is  a  good  approximation  for  jn,  retaining  only  the  resulting  Born 
term  in  (15)  as  the  first  order  approximate  solution,  and  ignoring  the  operation 
[1  -  (rn)]  ~  1  yields  the  following  result  for  the  total  normalized  current 

J„(l)  =  (V>  +  {*4,+T„-f„i}  (16) 

Reintroduction  of  the  phase  normalization  factor  yields  the  following  result  for 
what  has  been  called  the  normalized  first  order  smoothing  (NFOS)  approximation  [9]; 

j(’,1>  =  (J„i  +  r~.f„-p  (17) 

where 

<t>  =  exp  [Kk^-k^K]  (18) 

Taking  the  fluctuating  part  of  (17)  yields; 

*J{" 1 }  =  iJ1’  +  W  r\,  •  jV  +  64  r  (19) 

Comparing  (19)  with  (11),  i.e.  NFOS  with  FOS,  clearly  shows  that  (19)  contains  an 
additional  term  which  does  not  go  to  zero  as  k0^(C^)— >oo;  in  particular  this  term  is 

S+?n*ln  (20) 

If  it  turns  out  that  the  operator  involving  L(Pn)  in  (15)  must  be  retained  then  one  can 
deal  with  the  Fourier  transform  of  (15)  because  the  operator  becomes  a  simple 
multiplicative  factor  in  the  transform  domain.  Dealing  witii  the  transform  of  (15)  or 
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(19)  has  the  added  advantage  of  producing  a  quantity  that  is  proportional  to  the 
scattered  field.  In  the  transform  domain,  (20)  becomes 

#  *  (21) 

where  the  symbol  <g>  denotes  convolution.  This  term  has  an  important  interpretation 
which  will  be  discussed  later  in  the  paper. 

Capabilities  and  Limitations  of  Normalized  Smoothing 

The  limitation  of  the  FOS  approximation  is  that  it  will  fail  if  the  surface  height 
becomes  large  relative  to  a  wavelength.  If  the  surface  height  is  small,  FOS  has  no 
limitations  [8].  Comparisons  with  exact  numerical  solutions  will  be  required  to  establish 
the  height  range  over  which  FOS  exhibits  acceptable  accuracy.  However,  this  can  be 
accomplished  for  one  dimensionally  rough  surfaces  where  numerical  solutions  are 
available. 

Establishing  the  capabilities  and  limitations  of  NFOS  is  somewhat  more  difficult. 
To  accomplish  this,  the  following  approach  has  been  adopted.  The  NFOS 
approximation  is  based  on  the  neglect  o::  the  last  term  of  (15),  i.e. 

(22) 

The  only  time  this  term  is  known  to  be  of  negligible  consequence  is  when  the  surface 
height  is  small  because  this  is  when  NFOS  becomes  FOS.  Thus,  the  question  to  be 
addressed  is  when  does  the  term  in  (22)  look  like  it  corresponds  to  the  case  of  small 
height?  More  specifically,  the  question  should  be  when  does  the  propagator  <5rn  look  the 
same  as  with  the  small  height  FOS  approximation? 

The  height  dependence  in  6T,,  appears  as  the  difference 

AC  =  <(p)  -<0(po)  (23) 

both  in  the  phase  and  the  amplitude;  see  (4d),  (3d),  and  (3e).  Starting  off  with 

C  (p)  =  <8(p)  +  Z  (24) 
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where  Cs(p)  is  a  spatially  varying  small  height  component  and  Z  is  a  random  number 
yields 


t{p)  C  o(Po)  —  Cs(p)  Upo)  (2b) 

which  is  the  difference  of  two  small  surface  heights  and  inst  like  what  appears  in  FOS. 
The  decomposition  in  (24)  corresponds  to  the  sum  of  a  spatially  variable  small  height, 
<s(p),  superimposed  on  a  randomly  elevated  planar  surface,  Z.  Thus,  for  this  class  of 
surface  roughness,  NFOS  will  be  very  accurate.  Equation  (24)  suggests  that  the  next 
class  of  surface  to  be  considered  is  a  small  variable  roughness  superimposed  on  a 
randomly  elevated  and  tilted  planar  subsurface,  i.e. 

C(p)  =  Ca(p)  +  Z  +  m  •  p  (26) 

where 


m  =  mxx  +  myy 

and  mj  and  my  are  the  x  and  y  components  of  slope  of  the  planar  subsurface. 
Substituting  (26)  into  (23)  yields 

C(p)  —  Co(po)  =  ^s(^)  —  Cs(^o)  +  m  •  (^  —  p0 )  (27) 

The  propagator  <5r„  in  (22)  for  the  height  in  (27)  is  exactly  like  a  small  height  situation 
only  with  a  tilt  to  the  underlying  planar  surface.  However,  the  term  in  the  small  height 
case  corresponding  to  (22)  is  insensitive  to  global  tilting  of  the  underlying  surface.  This 
implies  that  (22)  is  negligible  for  the  case  of  small  scale  roughness  superimposed  on  a 
randomly  elevated  and  tilted  planar  subsurface.  It  should  be  noted  that  there  is  no 
restriction  on  the  magnitude  of  the  height  or  the  slope  of  the  subsurface. 

The  next  step  in  constructing  a  surface  to  which  NFOS  will  apply  is  to  add  .some 
underlying  curvature  to  (26),  i.e. 


<{p)  =  Cs(p)  +  Z  +  m»p  +  c»  R  (28) 
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where  c  is  a  vector  having  (constant)  curvature  components,  i.e. 


C  —  Cfx  X  +  Pyy  jy  Z 


and 


R  =  (x'2/2)  X  +  (y"/2)  y  +  (xy)  z 
The  height  difference  at  two  points,  using  (27),  will  be 

C (p)  -  Upo)  =  Up)  ~  <s(Po)  +  ifi  •  (p  -  Po)  +  c .  (R  -  R0)  (28) 

Unfortunately,  the  presence  of  the  curvature  terms  make  it  impossible  to  make  the 
NFOS  propagator  6r„  look  like  the  FOS  propagator  SV.  Hence,  when  the  underlying 
surface  has  a  nounegligiable  curvature,  it  is  no  longer  possible  to  rigorously  argue  that 
NFOS  is  accurate.  The  inability  of  NFOS  to  properly  account  for  a  nonnegligable 
underlying  surface  curvature  comprises  the  fundamental  limitation  of  NFOS.  However, 
it  should  be  noted  that  NFOS  provides  a  result  which  is  very  much  like  the 
conventional  composite  surface  scattering  model.  This  can  be  understood  by  examining 
(19).  The  First  term  on  the  right  hand  side  is  the  Kirchhoff  approximation  which 
dominates  about  the  specular  direction.  The  second  term  is  negligible  when  the 
underlying  height  is  large  because  (^)  =  0.  The  third  term  in  (19)  corresponds  to  the 
tilting  of  the  Rice  or  Bragg  scatter  solution  by  the  random  slopes  of  the  underlying 
surface.  Equation  (19)  is  an  improvement  over  the  composite  surface  model  in  that  the 
small  scale  surface  need  not  have  small  height  derivatives.  Conversely,  it  is  deficient  in 
that  it  permits  essentially  no  curvature  in  the  underlying  surface.  Thus,  it  can  be  said 
that  NFOS  provides  a  composite  surface  scattering  model  with  improved  accounting  for 
the  small  scale  surface  structure  but  no  large  scale  shadowing  effects.  Note  that  the 
primary  effect  of  the  large  scale  or  underlying  surface  structure  is  to  give  rise  to 
shadowing  in  the  conventional  composite  surface  model. 

One  possible  way  to  improve  on  (19)  is  to  iterate  (15)  one  or  more  additional 
times.  However,  it  is  not  clear  what  this  will  produce.  In  fact,  it  is  well  known  that 
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accounting  for  such  effects  as  shadowing  may  require  many  iterations  of  (15).  Hence, 
further  iteration  of  (15)  does  not  appear  to  be  a  profitable  effort.  In  the  next  section,  a 
somewhat  different  approach  will  be  developed  to  provide  an  improvement  to  the  NFOS 
approximation. 

A  Combination  of  The  Normalized  Smoothing  and  Composite  Surface  Ideas 

In  this  section,  ideas  suggested  by  both  the  NFOS  approximation  and  the 
composite  surface  model  will  be  combined  to  produce  a  more  robust  scattering  model. 
The  first  step  is  suggested  by  the  composite  model  and  entails  splitting  or  partitioning 
the  surface  roughness  spectrum  into  contiguous  parts  denoted  as  large  and  small  scale 
subspectra  as  shown  in  Figure  2.  The  partitioning  wavenumber  kp  is  chosen  to  satisfy 
two  criteria.  First,  it  must  be  less  than  the  electromagnetic  wavenumber  k0;  the 
smaller  it  is  relative  to  k0  the  better.  The  reason  for  this  is  to  insure  that  all  large  scale 
surface  height  spatial  frequencies  are  small  compared  to  k0.  This  condition  will  permit 
the  use  of  quasi-optical  approximations  in  calculating  the  scattering  from  the  large  scale 
part  of  the  surface.  It  is  implicitly  assumed  that  the  small  scale  spectrum  contributes 
most  significantly  to  the  total  surface  curvature,  rate  of  change  of  curvature,  etc.  It  is 
acceptable  for  the  large  scale  spectrum  to  contain  large  surface  slopes,  but  the  higher 
order  surface  height  derivatives  in  this  spectral  region  must  be  small.  The  second 
criterion  to  be  satisfied  by  the  choice  of  kp  is  that  the  FOS  approximation,  as  given  by 
(11)  with  C  =  Cs,  is  an  adequate  representation  for  the  scattering  from  the  small  scale 
structure.  This  implies  that  is  wot  too  large;  the  exact  range  is  yet  to  be 

determined. 

The  conditions  imposed  on  the  selection  of  kp  are  not  unlike  those  established  for 
the  composite  model  [12].  If  anything,  they  are.  probably  less  demanding  and  it  will  be 
subsequently  shown  how  they  may  be  relaxed  even  further.  In  any  case,  it  is  important 
to  note  that  many  natural  surfaces  do  exist  for  which  the  spectral  partitioning  may  be 
accommodated. 

The  next  step  is  suggested  by  the  NFOS  analysis;  equation  (2)  is  normalized  by 
the  phase  factor 

«*p[}(k** -*,-.)<:/]  (29) 

which  depends  only  on  the  large  scale  height.  Equation  (2)  becomes 
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Figure  2.  A  partitioning  of  the  surface  height  spectrum  into  large  and  small  scale 
parts.  Note  that  the  electromagnetic  wavenumber,  k0,  should  be  larger  than  the 
partitioning  wavenumber,  kp,  in  order  to  keep  the  large  scale  spectrum  smoothly 
varying  relative  to  X  0 . 
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(30) 


“  j»/‘  +  Lrn/  *  3«/ 

where 

Jn/  =  J  exp[-j(kS2  -  kiz) C/]  (31a) 

V  =  2N  *  H'(f)  expOk^C,  +  jkS2Cs]  (31b) 

rn/  =  ~  exp[jkS2A<s  +  jk^AQ]  N  x  [Ar  x  (31c) 

and  the  argument  of  the  free  space  Green’s  function,  as  a  reminder,  is 

[Ar2  +  Ay2  +  AC2] 1/2  =  [Ar2  +  Ay2  +  (ACS  +  AC/)2] 1/2  (32) 

When  dealing  with  the  NFOS  approximation,  it  was  found  that  NFOS  would  work  very 
nicely  for  an  underlying  surface  comprising  a  randomly  elevated  and  inclined  plane. 
The  way  to  adapt  this  result  to  (30)  is  as  follows.  The  integral  in  (30)  is  split  into  two 
contiguous  parts;  the  first  part  encompasses  some  yet  to  be  determined  neighborhood  of 
the  point  p0  =  p  while  the  second  part  includes  the  remainder  of  the  plane,  e.g. 

oo 

L  =  1 1  ( • )  dpo  =  1 1  ( • )  (ip0  +  |  ( • )  dp0  (33) 

~°°  \p-Po\  <Q  I P-Po  >Q 

Expanding  C/G  about  the  point  p0  =  p,  Q  is  determined  by  requiring  that  the  two  term 
approximation  for  C/0,  i-e. 

<10*  <1-  ^<lo\p0=p*^P  (34) 

is  accurate  when 

I &p\  -\p-Po\  <  Q  (35) 

Substituting  (34)  in  (32)  yields,  for  (35)  to  be  satisfied, 
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(36) 


I  At  |  =  [Ar2  +  Ay2  +  (A<s  +  C/jpAx  +  C/yAy)2]' /2 

Substituting  this  result  in  the  propagator  in  (31c)  shows  that  as  long  as  (35)  is  satisfied, 
the  integral  term  looks  like  the  case  of  small  scale  roughness  superimposed  on  a 
randomly  elevated  and  tilted  plane.  Hence,  the  integral  equation  in  (30)  can  be  written 
as 


where 


J 


Lr 

o 


(37) 


Jfl  =  V+Lri,,.J„,  (38) 

and 

L  =  ||  (*H?„  (39) 

\p-Po\  >Q 

L  =||  (•)<%  («) 

\p  ~  Po  I  •>  Q 

In  view  of  the  restricted  range  range  of  the  integral  in  (38)  and  the  fact  that  is 
constant  in  this  range,  it  seems  obvious  that  J ^  is  the  current  J  „/  with  the  large  scale 
slope  held  constant. 

Jfi  -  V  q  (41) 

However,  rather  than  approaching  the  solution  from  this  somewhat  heuristic  point  of 
view,  consider  the  following  approach.  According  to  (31b),  JfJj*  is  independent  of  ^ .  In 
addition,  by  virtue  of  the  way  the  area  I  in  (38)  was  constructed,  is  also  independent 
of  C/  in  (38).  Hence,  Jg  will  be  independent  of  the  large  scale  surface  height  .  This 
observation  implies  that  (37)  can  be  solved  by  iteration  in  which  each  integration  is 
performed  by  the  stationary  phase  approximation  provided  k,*2{C/2)  >  >  1.  Stationary 
phase  can  be  used  because  of  the  phase  factor 
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exp(jki2.AC,) 


in  rnt  (see  (31c))  and  because  the  point 

Po  =  P 

has  been  specifically  excluded  from  the  range  of  integration  in  (37).  The  solution  of 
(37),  under  these  conditions,  can  be  inferred  directly  from  previous  work  on  integral 
equations  of  this  kind  [13]; 

^nl  ~  ■-*/  3  b  +  roms  (42) 

where  S/  is  the  large  scale  dependent,  incident  shadowing  function  defined  to  be  unity  if 
a  point  on  the  surface  is  not  shadowed  from  the  incident  field  by  the  large  scale  surface 
and  zero  if  it  is.  It  is  important  to  note  that  only  the  large  scale  surface  enters  into  the 
determination  of  the  shadowing  function.  The  term  “roms”  in  (42)  stands  for  “ray  optic 
multiple  scattering”  and  it  represents  the  rays  that  bounce  around  on  the  surface  (by 
two  or  more  bounces)  before  traveling  into  the  upper  free  space.  Although  these  rays 
can  be  accounted  for,  it  is  a  most  tedious  “book  keeping”  effort  to  do  so  and  they  will 
be  ignored  in  the  remainder  of  this  paper,  e.g. 

3,i/  555  (43) 

Substituting  (43)  in  (37)  yields 

=  3n/  +  Ls/  rn/  •  JB  (44) 

It  is  interesting  to  consider  what  (43)  and  (44)  imply  relative  to  shadowing.  If  the  area 
I  is  shadowed  then  S/  is  zero  in  (44)  and  this  yields 

However,  substituting  this  result  in  (43)  produces 
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J„i = 0 


(45) 


Conversely,  if  tke  area  /  is  not  shadowed  so  that  S/  =  1  then  using  (43)  in  (44)  yields 

Jn/  =  V  +  Lrn/.Jfl/  (4G) 

Thus,  (45)  and  (46)  are  the  appropriate  solutions  when  either  the  point  (p,^  +  Cs)  is 
shadowed  or  illuminated,  respectively,  by  the  incident  field.  Consequently,  it  is  only 
necessary  to  multiply  the  solution  of  (46)  by  the  large  scale  shadowing  function  to  have 
a  complete  accounting  for  the  effects  of  shadowing.  Equation  (46)  can  be  solved  using 
the  method  of  smoothing  because  it  was  set  iq>  to  be  amenable  to  such  a  method.  The 
only  point  of  caution  is  to  recall  that  the  large  scale  surface  slope  is  taken  to  be 


constant  in  the  area  /.  Proceeding  as  in  obtaining  (16)  yields  1 

b((1)  =  V  +  p'„,.<J„/)  (47) 

or,  with  the  inclusion  of  the  shadowing  function, 

+  (48) 

where,  in  review, 

J„/1  =  2ri0N  x  hj  exp(j(k„  -  kizKs  -  jka  •  p]  (49a) 

=  2H0  (N|  V  q  X  hj-  (exp[j(kS2  -  k}>)Cs])exp(  -  jkit  •  p)  (49b) 
r„/  =  -  exp(jkSiACs  +  jklz  V  C,  •  Ap]  N  X  [Ar  x  (49c) 


’It  is  assumed  that  the  Fourier  transform  of  the  average  propagator  (T  nl)  is 
negligibly  small.  If  this  is  the  case  then  L^rn^  may  also  be  ignored.  If  the  average 
propagator  cannot  be  ignored,  it  can  be  shown  that  (47)  becomes 

•4/(1)  =  (>  -  L<r„,>]  - 1  (4/  +  oiy,  -  L(r„,)j  -  ■<.)„/» 
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ACs  —  (s  (so 


(49d) 

Ap  =  (x-x0)x  +  (y-y0)y  =  Axx  +  Ayy  (49e) 

V  C/  =  d^/dx  x  +  (Ktfd y  y  (49f) 

Ar  =  {Ax2  +  Ay2  +  {ACS  +  V</- A£)2} 1 I2  (49g) 

Also,  H0  is  the  amplitude  of  the  incident  magnetic  field,  li,-  is  its  direction,  and  ^N|  y  ^ 
is  the  conditional  average  of  N  conditioned  on  holding  the  large  scale  slopes  constant.  It 
must  also  be  remembered  that  in  (48),  the  large  scale  slopes  V(,  are  to  be  held 
constant.  To  complete  the  evaluation  of  the  current,  the  normalizing  phase  factor  must 
be  reintroduced,  e.g. 


J  -  s/  </>  { +  Lrn/  •  (50) 

where 

<j>  =  exp(j(ksz  -  kjZXj)  (51) 

The  far  zone  electric  field  scattered  by  the  surface  is  given  by 

ES(R)  =  -iK'loSW  {?  (£„)-[£,.?  (ks,))ks)  (52) 

where  R  is  the  distance  from  the  origin  of  the  coordinate  system  on  the  mean  surface  to 
the  point  of  observation,  ks  is  a  unit  vector  pointing  in  this  direction, 


R  =  Rks 

and 

J  (Kt)  =  |  J(*()  exp(j  kst  •  p)dp  (53) 

where 

ksi  =  ko(sin0s  cos  <t>s  x  +  sin 0S  sin  <f>s y) 
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p  =  XX  +  yy 


Discussion  of  Results 

Equation  (50)  is  the  central  result  of  this  paper.  To  the  author’s  knowledge,  it 
represents  the  first  tractable  description  of  the  current  induced  on  a  multiscale  surface. 
There  are  certain  advantages  to  dealing  with  the  current  rather  than  the.  far  zone 
scattered  field.  This  becomes  especially  clear  if  it  is  desirable  to  study  Fresnel  zone 
scattering  from  a  rough  surface. 


While  (50)  contains  a  number  of  approximations,  it  represents  a  considerable 
improvement  over  the  conventional  composite  surface  scattering  model.  Paramount 
among  these  is  the  fact  that  the  small  scale  part  of  the  surface  need  not  have  small 
surface  height  derivatives.  This  result  follows  from  the  use  of  FOS  to  estimate  the 
current  induced  on  the  surface  within  the  surface  area  I.  In  the  conventional  composite 
model,  boundary  perturbation  theory  is  used  to  estimate  the  effects  of  the  small  scale 
structure  and  it  is  much  more  limited  in  its  range  of  applicability.  Another  subtle 
improvement  provided  by  (50)  is  use  of  a  finite  area,  /,  rather  than  an  infinite  one. 
This  finite  area  is  tailored  to  the  large  scale  slopes.  That  is,  it  is  the  largest  area  over 
which  the  effects  of  large  scale  curvature  may  be  ignored.  Thus,  if  one  were  dealing 
with  surfaces  for  which  the  large  scale  rms  slopes  are  not  too  large  but  large  scale  rms 
curvature  is  moderate2  then  I  in  (50)  would  have  to  be  decreased  relative  to  its  value 
when  the  rms  slope- to-curvature  ratio  were  larger.  It  should  be  noted  that  this  effect 
builds  into  the  solution  a  dependence  upon  the  large  scale  curvature  that  is  quite 
independent  of  the  method  used  to  analyze  the  scattering  from  the  large  scale  surface 
structure.  This  dependence  requires  further  study  in  order  to  understand  its  full 
importance.  However,  it  is  possible  to  provide  some  estimate  of  the  effect  of  reducing  / 
as  follows.  If  /  is  infinite, 


r„/(£-kii)  •]*„/* 


(54) 


2It  is  assumed  that  the  large  scale  curvature  is  sufficiently  small  to  permit  the 
use  of  stationary  phase  integration  techniques  on  the  large  scale  surface. 
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where  the  tilde  denotes  the  Fourier  transform  (from  pQ  to  -£(<)  and  j  ^  is  a  constant 
vector.  If  I  is  finite,  it  is  necessai;,  to  convolve  (54)  with  a  transformed  support 
function  representing  the  finite  size  of  /.  This  convolution  has-  the  effect  of  blurring  or 
smearing  the  sharply  defined  effects  contained  in 

It  is  interesting  to  see  how  (50)  reduces  to  the  conventional  composite  surface 
scattering  model.  The  details  are  very  involved  and  will  only  be  highlighted  here.  The 
term 


in  (50)  is  what  gives  rise  to  the  quasi-optical  part  of  the  scattering  that  is  dominant  in 
and  about  the  specular  direction.  Note  that  this  term  contains  large  scale  incidence 
shadowing  (S/),  deep  phase  modulation  associated  with  the  large  scale  height  (see 
equation  (51)),  and  attenuation  due  to  the  small  scale  height  ((•I,,/))-  Apart  from  this 
small  scale  dependent  attenuation,  the  remaining  factors  are  large  scale  dependent  and 
have  been  treated  previously  by  Saucer  [14].  If  the  small  scale  surface  height  is  very 
small  and  the  large  scale  curvature  is  so  small  that 


then  the  terms 


SJV  +  Jr,  i, 

in  (50)  will  go  to  the  Rice  boundary  perturbation  solution  on  a  tilted  plane;  the  tilt 
being  provided  by  the  large  scale  surface  slopes.  After  multiplying  by  S/d,  taking  the 
Fourier  transform  to  generate  the  scattered  field,  squaring  and  averaging  to  produce  the 
second  moment  of  the  scattered  field,  the  net  effect  is  to  obtain  a  smeared  Bragg 
scattering.  The  smearing  results  from  the  range  of  possible  values  for  the  large  scale 
slopes  as  dictated  by  the  probability  density  function  of  the  slopes.  Thus,  (50) 
reproduces  the  conventional  composite  surface  model  under  the  same  set  of 
assumptions  that  are  usually  associated  with  the  conventional  model. 

There  is  one  final  point  about  the  method  introduced  here  that  warrants 
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discussion.  While  it  is  clear  that  a  small  scale  surface  can  always  be  split  off  from  the 
total  surface,  it  is  not  so  obvious  that  the  resulting  large  scale  features  will  satisfy  the 
criteria  placed  on  them.  For  example,  it  may  be  necessary  to  include  wave  diffraction 
effects  in  the  solution  of  (37)  rather  than  just  wave  reflection  and  shadowing  as  in  (43). 
In  a  purely  formal  sense,  the  solution  of  (37)  is  given  by 


Ui/  —  [•  -  Lru/]  1  •  J  w 


so  the  integral  equation  for  .J  ^  is 


(55) 


(56) 


and  within  the  confines  of  first  order  smoothing 


*JB(1)  =  (l-(A)r‘.  {U1(f  +  OA  •  [l  -  (A)]"1  •<•?„/)}  (57) 


where 


A  = '%./  *  I1  -  ^r„,;  - 1 

Thus,  (55)  becomes 

•b/1’  =  [I  -LI'J-  1  .  {[I  -(A)]  -  1  .(.U)  +  «JB(1)}  (58) 

Equation  (58)  represents  an  approximate  solution  for  the  normalized  current  based  on 
an  exact  solution  for  the  long  range  interactions  on  the  surface  and  a  first  order 
smoothing  approximation  for  the  short  range  interactions.  Equation  (55)  represents  the 
former  while  (57)  represents  the  latter.  It,  is  clear  from  (55)  and  (58)  that  the  long 
range  solution  is  highly  dependent  upon  ones  ability  to  invert  the  operator  [1  -Lr^].  In 
the  optical  limit,  (43)  shows  that 

(5!» 

where  ray  optic  multiple  scattering  terms  have  been  ignored.  If  wave  diffraction  is  to 
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be  included,  it  will  be  necessary  to  do  something  like  suggested  by  Chaloupka  and 
Meckelburg  [15]  and  Ansorge  [16],  The  important  point  is  that  the  new  method 
presented  in  this  paper  not  only  improves  on  the  classical  composite  surface  scattering 
model  but  it  clearly  shows  how  further  improvements  must  proceed. 

Summary 

This  paper  adapts  and  combines  the  fundamental  principles  comprising  the 
conventional  composite  surface  scattering  model  and  the  normalised  first  order 
smoothing  (NFOS)  approximation  to  provide  a  new  approach  to  the  problem  of 
scattering  by  multiscale  surfaces.  The  conventional  composite  surface  model  suggests 
partitioning  the  surface  into  two  subscales  which  have  rms  heights  that  are  small  and 
large  relative  to  a  wavelength.  The  NFOS  approximation  suggests  normalizing  the 
basic  integral  for  the  current  induced  on  the  surface  by  the  large  scale  height  dependent 
phase  factor  contained  in  the  Kirchhoff  term  to  produce  an  integral  equation  that  is 
more  amenable  to  the  first,  order  smoothing  approximation. 

The  real  heart  of  the  improvement  rendered  by  this  approach  stems  from 
splitting  the  integral  in  the  current  integral  equation  into  the  sum  of  one  encompassing 
the  immediate  neighborhood  of  the  observation  point  and  one  including  all  of  the 
remaining  surface.  The  extent  of  the  former  integration  is  chosen  such  that  the  large 
scale  surface  has  essentially  a  constant  slope  within  this  region.  The  sum  of  this 
integral  and  the  Kirchhoff  term  act  as  the  Born  term  in  the  integral  equation  involving 
the  remaning  integral  term.  This  equation  may  be  solved  by  iteration  in  the  high 
frequency  limit  and,  hence,  yields  the  shadowed  Born  term  as  the  solution.  Thus,  the 
net  effect  of  the  large  scale  solution  at  this  point  is  simply  to  shadow  the  new  Born 
term.  The  solution  of  the  resulting  finite  range  integral  equation  is  accomplished  by 
first  order  smoothing  and  is  valid  as  long  as  the  small  scale  height  is  not  too  large 
(relative  to  the  wavelength)  but  its  validity  is  independent  of  lie  surface  slopes, 
curvature,  etc. 

The  net  result  of  this  analysis  is  an  expression  for  the  surface  current  which 
contains  all  of  the  fundamental  features  of  the  conventional  composite  surface  scattering 
model  plus  some  marked  improvements.  The  latter  entail  accounting  for  more  general 
small  scale  surface  structure  and  an  estimate  of  the  effects  of  the  large  scale  surface 
curvature  on  the  scattering  from  the  small  scale  features.  Finally,  this  result  is 
amenable  to  correction  for  such  phenomenon  as  high  frequency  diffraction  by  the  large 
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scale  surface.  Although  it  is  essential  to  expose  this  result  to  extensive  computation  for 
the  purpose  of  establishing  its  range  of  validity,  the  basic  merits  of  the  new  result  are 
clear.  Of  particular  note  is  the  fact  we  now  have  an  expression  for  the  current  induced 
on  a  multiscale  rough  surface  rather  than  just  one  of  scattered  field  moments.  The 
advantages  of  this  knowledge  are*  obvious. 
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